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Abstract 



We study the time and temperature dependent correlation functions for an impenetra- 
ble Bose gas with Neumann or Dirichlet boundary conditions {'il}{xi,0)ip^ {x2,t))±^T- 
' We derive the Fredholm determinant formulae for the correlation functions, by means 

Q> \ of the Bethe Ansatz. For the special case xi — 0, we express correlation functions with 

0^ ' Neumann boundary conditions ('(/'(O, O)'0^(x2, t))j^^Ti in terms of solutions of nonlinear 

(— I \ partial differential equations which were introduced in as a generalization of the 

^ I' nonlinear Schrodinger equations. We generalize the Fredholm minor determinant for- 

mulae of ground state correlation functions (V'(2;i)'(/'^(^2))±,o in 0? to the Fredholm 
determinant formulae for the time and temperature dependent correlation functions 
(V'(a;i,0)V't(a:2,t))±,T, <GR, 'r>0. 



1 Introduction 



> 

X 

' In the standard treatment of quantum integrable models, one starts with a finite box 

and impose periodic boundary conditions, in order to ensure integrability. Recently, there 
has been increasing interest in exploring other possible boundary conditions compatible 
with integrability. These other possible boundary conditions are called "open boundary 
conditions" . 

With open boundary conditions, the works on the two dimensional Ising model are 
among the earliest. By the help of graph theoretical approach, B.M. McCoy and T.T. Wu 
1^] studied the two dimensional Ising model with open boundary conditions. They cal- 
culated the local magnetizations. E.K. Sklyanin began the Bethe Ansatz approach to 
open boundary problems. M. Jimbo et.al. ||5| studied the antiferromagnetic XXZ chains 
with open boundary conditions and derived an integrable representation of correlation 
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functions, using Sklyanin's algebraic Bethe Ansatz framework and representation theoret- 
ical approach invented by Kyoto school [P, 0]. T. Kojima ^ studied the ground state 
correlation functions for an impenetrable Bose gas with open boundary conditions: 

{^p{xl)^p\x2)). 

Kojima derived the Fredholm minor determinant representations for the ground state 
correlation functions by the help of fermions, which have the integral kernel: 

sin(A ~ /^) _|_ sin(A + fi) 
X — fi A + /i 

The integral intervals depend on the space parameter xi, X2- In this paper we study an 
impenetrable Bose gas with open boundary conditions. We are interested in the finite- 
temperature dynamical correlation functions: 



{^{xi,QW{x2,t)) 



We derive the Fredholm determinant representations for the dynamical correlation func- 
tions by the coordinate Bethe Ansatz, which have the integral kernel: 

(i(A, ^) ± L(A, -m)) 

where we have used 

L(A,^) = ^- ie**^ sin(xi(A - //)) + e**^ sin(x2(A - /x)) 

(1.1) 

+1 P.V. r ( — e^*"' sin((s - /i)xi) sin((s - X)x2)ds\ . 

71" J-oc\s-f^ s-XJ J 

Here the notation P.V. represents Cauchy's principle value and the measure '!9(A) is 
given by 

m = W-v- (1-2) 



1 -|- exp 



X2-h 

T 



For an impenetrable Bose gas without boundaries, V. Korepin and N. Slavnov |^] has 
derived the Fredholm determinant formulae for the dynamical correlation functions. 

To describe the 2n point dynamical correlation functions for an impenetrable Bose gas 
without boundaries, N. Slavnov [Q] introduced a system of nonlinear partial differential 
equations, which becomes the nonlinear Schrodinger equation in the simplest case. The 
generalization of the nonlinear Schrodinger equations has 2n time variables tj, (1 < j < 2n) 
and 2n space variables xj, (1 < j < 2n). In this paper, we consider the dynamical 
correlation functions with Neumann boundary conditions for the special case that one 
space parameter has the value xi = 0: (^(0, 0)^/^^(x, We express the dynamical 

correlation functions in terms of a solution of Slavnov's generalization of the nonlinear 
Schrodinger equations. The differential equations, which describe four-point correlation 
functions without boundaries: 



{lp{xi,ti)'lp\x2,t2)lp{x3,t3)lp\x4„ ti)) 
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describe the dynamical correlation functions with Neumann boundary conditions: 



Now a few words about the organization of the paper. In Section 2 we formulate the 
problem and summarize the main results. In Section 3 we obtain the determinant formulae 
for the field form factors. In Section 4 we obtain the Fredholm determinant representation 
for the dynamical correlation functions. In Section 5 we consider the completely integrable 
differential equation which describes the 2-point dynamical correlation functions with 
Neumann boundary conditions. In Section 6 we consider the special case that time t = 
and derive the Fredholm minor determinant representations for the finite-temperature 
fields correlation functions. We show that our Fredholm formulae coincides with the one 
which has been obtained ||2| at temperature T = 0. 

2 Formulation and Results 

The purpose of this section is to formulate the problem and summarize the main results. 
The Hamiltonian of our model is given by 



and L > is the size of box. The parameters h > and /iq S R represent the chemical 
potential and the boundary chemical potential respectively. We only consider the case 
of the coupling constant c = oo, so-called "impenetrable case". The Hamiltonian H acts 
on the Fock space of the Bose fields defined by the following relations between the Fock 
vacuum |0) and the Bose fields: 



where the integrand ip^i^i, . . . , zj\f) is a C-valued function. The eigenvector problem: 
H\^]\f) = E]\f\^]y), {E]\f G R), is equivalent to the quantum mechanics problem defined 
by the following four conditions of the integrand function ■ip]\f{zi, . . . , zn)- 

1. The wave function ipN = ipN^zi, . . . , zn) satisfies the free-particle Schrodinger equa- 
tion in the case of variables < / Zj < L: 



(V'(o,o)VHa:,t)>T,+. 




(O|'0"^(x) = 0, ij{x)\0) = 0, (0|0) = 1. 
A A^-particle state vector I^'at) is given by 





tPn{zi, . . . ,zn) = En ■ TpNizi, ...,zn) 



{0 < Ziy^ Zj < L, En e R). 
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2. The wave function is symmetric with respect to the variables: 

3. The wave function ip^ satisfies the integrable open boundary conditions: 



dz.j 



0, 

:0, (j = l,...,A^). 



4. The wave function ip^q vanishes whenever the coordinates coincide: 

iI}n{zi, ...,Zi,...,Zj,..., ZAr)|^_^^. = 0. 
This condition corresponds to the condition: c — > oo. 

The wave functions tpN which satisfy the above four conditions were constructed 
They are parameterized by the spectral parameters 

iPn{zi,...,zn\Xi,---,Xn) 

= Cons. T 1 sgn(2;j — zj.) det {Xj cos(AjZfc) + ho sm{XjZk)) ■ 

l<j,k<N 
i<j<k<N - 

X 

Here the function sgnfa;) = - — r and the spectral parameters < Ai < A2 < • • • < Xn are 

\x\ 

determined by the so-called Bethe Ansatz equations: 

Xj = jlj, {IjeN, j = l,2,...,N). (2.1) 

Because the coupling constant c ^ 00, the Bethe Ansatz equations become simple. The 
constant factor ^^Cons." is determined by 

(^^(Ai, . . . , A^)|^'7v(Ai, . . . , Xn)) = {2L f. 
The eigenvalue Ej\f{{X}): 

H\^n{Xi, . . . , A^)), = ^^({A})|^^(Ai, . . . , XN))e, 
is given by 

TV 

EN{{X})=Y.i^]-h). 
i=i 

We assume that the set {|^Af (Ai, . . . , AAr))}aii{A}jv A^eN is a basis of physical space of 
this model. Here the index all{A}7v represents all the solutions of the Bethe Ansatz 
equations ( p. 11 ). This type assumption is usually called "Bethe Ansatz". The following 
lemma is a foundation of our analysis. 
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Lemma 2.1 If the boundary condition /iq takes the special value ho = 0,oo, the eigenvec- 
tors |\I'({A})) satisfy orthogonality relations 



N 



(^Ar(Ai,...,Ajv)|^'w(/ii,...,m)) = (2L)^ J]5a,„ 

i=i 

Here 5x,^ is Kronecker Delta. 



{ho = 0,00). 



(2.2) 



To prove the above lemma, we have used the Bethe-Ansatz equations of the spectral 
parameters. In the sequel we use the orthogonality relations of the eigenstates, therefore 
we concentrate our attentions to the case of the special boundary conditions: ho = 0, 00. 
The boundary conditions ho = and ho = 00 are called Neumann, Dirichlet, respectively. 
In the sequel we use the following abberivations 



|'I'Ar(Ai, . . . , AAr))4_ for Neumann, 
The constant ^^Cons." is given by 

' N 



|^'7v(Ai 



,Xj\f))- for Dirichlet. 



Cons. 



1 



2i\ 



N\ \hoJ 



N 



for Neumann, 



for Dirichlet. 



By using the orthogonal relations ( |2.2| ) and the so-called "Bethe Ansatz" , we arrive at 
the completeness relation: 



id 



|^jv(Ai, . . . , Ajv))e e(^jv(Al, . . . , Ajv)| 

.(^-TvlAi, . . . , Aiv)|^iv(Ai, . . . , Ajv))e ■ 



(2.3) 



N=0 all{A}jv 

The Bose fields 'ip{x,t), tp^Xjt) are developed by the time t by 

idti^ = [i;,H], idti^^ = [i;\H]. 
More explicitly the time dependence of the Bose fields are written by 



ip{x,t) 



JHt 



'4j{x)e 



-iHt 



i;^x,t) 



(x)e" 



In this paper we are interested in the dynamical correlation functions ti)V'^(x2, t2))e,T 
defined by the following way. For the nonzero temperature T > 0, the dynamical correla- 
tion functions for state are defined by the summation of the every states: 



{'>p{xi,ti)lp\x2,t2))e,N,T 

En{{\})\ 




[all{A}jv 



T ) 

e(^jv(Ai, . . . , Ajv)|^(xi,ti)V;t(3;2, t2)|^jv(Ai, . . . , Ajy)) 
e(^'7v(Ai, . . . , A7v)|^'7v(Ai, . . . , Atv)). 



i^7v({A}) ^ 
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where the index e = it represents the boundary conditions. Now the index "+" represents 
the dynamical correlation functions with Neumann boundary conditions. The index "— " 
represents the dynamical correlation functions with Dirichlet boundary conditions. For 
the ground state case T = 0, the dynamical correlation functions for N state are defined 
by the vacuum expectation value of the ground state: 

{lp{xi,ti)'lp\x2,t2))e,N,0 

e(*7v(Al, . . . , XN)\i^{xi,ti)llj^X2,t2)\'ifN{Xl, • • • , A^)), 



e(^'^(Ai, . . . , Ajv)|^'iv(Al, . . . , A7v))e 

where the spectral parameters (Ai, . . . , Aat) are given by 



— (j — 1), for Neumann, 
Ij 

vr 

— j, for Dirichlet. 
L 



In this paper we are interested in the thermodynamic limit of the correlation functions, 
for the temperature T > and the time ti,t2 G R- For the nonzero temperature T > 0, 
the dynamical correlation functions in the thermodynamic limit are defined by 

{'tp{xi,ti)'>p\x2,t2))e,T = lim {'lpixi,ti)tp^{x2,t2))e,N,T- 



f =_D{T) 



Here the density D{T) = — Q is given by 

D{T) = - f ^{X)dX, (2.4) 
Jo 

where the Fermi weight ^{X) is defined in ( |1.2| ). For the ground state T = 0, the dynamical 
correlation function in the thermodynamic limit is defined by 

{lp{xi,ti)'lp\x2,t2))e,0 = lim {lp{xi,ti)lp\x2,t2))e,Nfi- 

N,L—foo 
f-=-D{0) 

N 

Here the density D[0) = — can be chosen arbitrary. In this paper we give the Fredholm de- 

Ij 

terminant representations for the dynamical correlation functions {^|J{xl,tl)^lJ^ {x2,t2))e,T, 
(T > 0, e = lb). Because the following relation hols: 

(V'(xi,0)'(/'^(x2,t2 - tl))e,T = {'4'{xi,ti)llj\x2,t2))e,T, 

we only need one time parameter t = t2 — ti, to describe correlation functions. In the 
sequel, we use the abbreviation t = t2 — ti. Let us set 

T{s\x,t) = its'^ — ixs, (2-5) 

G{x) = — / e"(^l^'*)ds, (2.6) 

and 

P(X\xi,X2) = /e^(^l^i'*) - - P.V. H ^—e^^^'^i'*) sin(x2(s - X))ds\ .(2.7) 

I vr J-ooS-X J 

In Section 4 we derive the following formulae. 
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Theorem 2.2 In the thermodynamic limit N,L — > oo, such that 'J^ ~ ground 
state dynamical correlation functions have the Fredholm determinant representations 

{i>{xi,0)lP\x2,t)),,0 

= e-*^* (g{xi - X2) + eG{xi + X2) + det (^1 - ^V, - aA,^ 

where the function G{x) is given in ( \2.(\ ). Here the integral operators and are 
defined by 

{VJ){X)= rK(A,/i)/(/i)d/i, {AJ){X)= /%,(A,^)/(^)d/i, 
JO Jo 

where the Fermi sphere q = ttD and the integral kernel are given by Neumann or Dirichlet 
sum: 

K(A,//) = L(A,^)+eL(A,-^), 

A^{X,fj.) = e(P(A|xi,a;2) + eP{-X\xi, X2)) {P{n\x2,xi) + eP{-fi\x2,xi)) . 

Here we have used the function L{X, fi) defined in and the function P(A|xi, X2) defined 
in ( \2. Here we can choose the density D > arbitrary. 



We have succeeded to write the integral kernel by elementary functions: 

roo „T{s\y,t) 

P.V. / -ds, 

J -00 S — A 

and trigonometric functions. In Section 4, we consider the finite temperature case, too. 

N 

Theorem 2.3 In the thermodynamic limit: N,L ^ 00, such that — = D(T) i\2.4) , the 
finite temperature dynamical correlation functions have the Fredholm determinant repre- 
sentations 

(V(xi,0)V^t(^2,t)).,r = e-'""' (g{xi - X2) + eG{xi + X2) + 
X det f 1 - -K T - T I 
Here the temperature T > and the integral operators V^^t A^^t o,re defined by 

roo /"OO 

{%Tf){x)= v,{x,f,)^{i^)f{i^)df,, (1,,t/)(A)=/ A,{x,^l)^^I)f{^^)d^I, 

Jo Jo 



where the Fermi weight '&{X) is given in fuTaj. 
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In Section 5, we derive the differential equations for correlation functions for the case 
xi = 0. To describe 2n point dynamical correlation functions for an impenetrabel Bose gas 
without boundaries, N. Slavnov introduced a system of nonlinear partial differential 
equations, which becomes the nonlinear Schrodinger equation in the simplest case. In 
this paper, we express the dynamical correlation functions {^p{0,0)'^p^ {x,t))T,+ , {T > 0) in 
terms of Slavnov's generalization of the nonlinear Schrodinger equations For xi = 
and Dirichlet boundary case: 

(V(o,o)V't(^2,t))T,- = 0, 

because the wave functions become to zero. We consider the case xi = and Neumann 
boundary conditions. First we consider T = case. Let us set 

(wf) (A) = w{x, ^I)f{^i)d^i, {q = ttD), (2.8) 

^ sm(x(A-/.)) ^ sin(x(A + ^))_ 
A — /i A + /X 



Theorem 2.4 The correlation functions for an impenetrable Bose gas with Neumann 
boundaries at the ground state are given by the following formulaes: 

, / —X X 
t t 



(V'(0,0)V^(x,t))o,+ = 2e-^^*det (^1 - ^W^ 61,4 ( 



Here the integral operator W is defined in (2.1 ) and the function 61^4 is a component of 
matrix b defined in (5.5). 

The matrix b defined in ( ^.5D satisfies a set of partial differential equations introduced 
in 0: 



—Lk - —M, + [Lk, M,] = 0, (1 < j, k<i). (2.10) 
otj onk 



d ^ d 
otj 

Here we have used 

db 

Lj ifi) = fiPj + [b, Pj], M, ifi) = -fiL, (/^) + ^ ' 
where we have used the matrix Pj whose components are defined by 

iPj)l,m = i^l,jSm,j- (2.11) 

The diffential equations ( p. 10 ) describe the logarithmic derivatives of the four point cor- 
relation functions without boundaries, too: 

In Section 5, we consider the finite temperature case, too. Let us consider the finite 
temperature case T > 0. Let us set 

(Wrf) (A) = Wt{X, fi)f{fi)dfi, (2.12) 

WT{X,fi) = WiX,fiMf,), (2.13) 
where the kernel W{X,^) is defined in (p.9[) and ^{^) is defined in (|1.2|). 
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Theorem 2.5 The correlation functions for an impenetrable Bose gas with Neumann 
boundaries are given by the following formulaes: 



(V'(o,o)vna:,t)> 



2e 



-iht 



det ( 1 - -Wt ] 6^4 



TT 



—X X 
t t 



Here the integral operator Wt is defined in (2.1^) and the function bj^ is a component of 
the matrix b'^ defined in 



The matrix 6^ satisfy a set of diffential equations ( |2.10| ), too. (We substitute b to b'^ .) 

T. Kojima [|| derived the Fredholm minor determinants formulae for the ground state 
correlation functions: (■0(xi)'0^(x2))o,e- In Section 6 of this paper, we consider the special 
case for the time t = of our Fredholm determinant formulae: 

(V(xi,0)V^+(x2,0))t,., (e = ±) 

and derive the Fredholm minor determinant formulae for temperature T > 0. This Fred- 
holm minor determinant formulae for T = coincide with the one which has been ob- 
tained Q. Let us set 



(eiy^f) (e) = ((E(2/i - e) + E(2/2 - e))Oe,T{^, ^o) fi^d^', (2.14) 



where 



■d{i') {cos((^ — r])v) + e cos{{£, + r])v)} du, 



and i?(A) is defined in ( |l.2D . Here E(^) represents the step function 

1, for ^ > 0, 
0, for C < 0. 



Theorem 2.6 For the temperature T > 0, the field correlation functions have the first 
Fredholm minor determinants representations 



1 



X2 
Xi 



(V(xi)^t(^2))e,T = ^det( I- -I 
where the integral operator ^^'^^^ is defined in (2.14)- 



Theorem 2.7 ^ The ground state correlation functions have the first Fredholm minor 
determinant formulae 



{^p{xi)^l)\x2))e,0 

Here the integral operator is defined by 



-det ( 1 - 
2 V vr ' 



X2 
Xl 



where 



N 



Here the density ^ ~ ^ chosen arbitrary. 
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3 Form Factors 

The purpose of this section is to derive the determinant formulae for the form factors. 
First we prepare a lemma. 

Lemma 3.1 For the sequences {fj^k}j=i,...,N+i, k=i,...,N oind {5(j}j=i^...^7v+i) the following 
holds: 



^ ( d\ 



.(3.1) 

a=0 



Proof. Consider the coset decomposition: 

Sn+1 = Sn{N + 1) U SNiN) ■ {N,N + 1)U---U SNil) ■{l,N + 1), 

where SnU) is permutations of (1, . . . , j — 1, + 1, j + 1, . . . , N). Rewrite the left side of 
the equation with respect to the coset decomposition: 

N N 
(L.H.S.) = ^ 5i X] '^Snlr • (j, iV + 1)) JJ fr{k),k ■ /r(Af+i) j 

N 

+gN+i Yl ^gn{T-{N + l,N + l))l[U(^k),k = {R-H.S.) 

T€SNiN+l) k=l 

Q.E.D. 

Now let us consider the field form factor: 

,(*Ar+l(Ai, . . . , AAr+i)|V'^(2;)|^'Ar(;Ui, . . . , ^Ar))e 

= VA^ + 1 / dzi... dZNll^*N+li^l^---^ZN,x\Xl,...,XN+l) 

Jo Jo 

= ^ y sgn (j(e"^^-(^+i)^ + ee^^-e^+i)^) 

V(1 + 5m,o)(1 + ^..o) 

]J I / sgn - x){e-'^-(i)^ + ee*^-(^)^)(e^'^^^ + ee-*''^'^)| . 
,=1 '-•^0 J 



N 

X 



To derive the third line, we have used a simple fact: 

N N 

Yl (TT Yl U(j),rij) = iV! Sgn a JJ f^^j^j. 

o-jTeSiv j=i o-e5]v i=i 



Using lemma 3.1, we arrive at the determinant formulae for the form factors. 
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Lemma 3.2 The field form factors have the determinant formula 

( d 



X det {Ie{x\Xj,fik) - aC^{x\\j)Ie{x\XN+i,fJ.k)) 

l<j,k<N 



a=0 



Here we have used 



C,{x\X) = ^ ^ (e-^^^ + ee^^^), (3.2) 



1 r 4 

^(1 + <5a,o)(1 + '5/.,o) U-M 

4 1 

+ e— - — sm(x(A + ^)) - 2L(5a,m + e^A,o'5/x,o) } ■ 



(3.3) 



We can write the field form factors without using integrals. 

From the relation -0^(3;, t) = e^^^ijj^ {x)e~'^^^ , the dynamical form factors are given by 

,(^'Ar+l(Ai, . . . ,XN+l)\'tp'^{x,t)\^N{l^l, ■ ■ • ,/"iv))e 

= exp|.t||-/i + ^'A|_^^2j| (34) 

X ,(^'Ar+i(Ai, . . . , Xn+iM\x)\^ n{1^1, ■ ■ ■,fJ'N))e- 

4 Correlation Functions 

The purpose of this section is to derive the Fredholm determinant formulas of the dynami- 
cal correlation functions {■ip{xi,ti)'ip^ {x2, i2))e,r- First we consider the vacuum expectation 
values of fields operators. Using the completeness relation (|2.3D, the vacuum expectation 
values of two fields are given by 

e(^'Ar(/il, . . . ,flN)\lPixi,ti)lp^{x2,t2)\'^NifJ'l, ■ ■ • ,/^Ar))e 



E 



(^7V(/^1, • • • ,m)I^Af(Atl>- • ■,m))e 

(M/^(M)|V(xi,ti)|*^+l({A}))e e(*iV+l({A})|V^(x2,t2)|*7v(M)>. 



,(M/;v(M)|^iv(M)>ee(^7V+l({A})|W({A}))e 
all{A}jv+i 



Using the equation ( |3.4| ) and the following relations: 

,{^Nmmxi,ti)\^N+ii{x}))e=e {^N+ii{xmHxi,ti)\^Nm)): 

,(^^({A})|^^({A})), = {2Lf, 
we obtain 



(iV + l)! \2LJ 



AiGfN Ajv+iefN 



X ,(*^+i({A})|^T(^^)|,i;^({^}))* ,(.I;^+i({A})|V^T(^2)|^^({M})). 
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The translation invariance of time holds 

,(^'iv(M)|V'(xi,tl)V^(x2,t2)|1'jv(M)). 

In the sequel we set the abberiviation t = t2 — ti. Remember a following simple fact. For 
sequences {/ji,...j„}ji,...j„e7, {5ji,...j„}ji,...j„G/, (/: some index set), the following holds 

(Sym /)j,„„j„(Sym5)j,,...j„ = ^ /ji,...,i„ (Sym 

Here we have used 
(Sym f)j 

CT(l)>---J(T(n) • 



The form factors have the determinant formulae in lemma 3.2 and 

AT 

o-€Sn+1 i = i 

We obtain 



JV 



2L J 

X det ( e**'*'j Ie(j;i|Aj,/Xfc)/e(j;2|Aj,/Xj) — ae^^'^i C*{xi\\j)If^{x2\Xj, Hj) 

l<j,k<N \ 



l<j,fc<Ar 

X e^*^^+i Ce (X2 1 Atv+I )/e (^1 1 A^V+l , /ifc ) 



a=0 



The J th line of the above matrix only depends on \j not on Xk,{k ^ j), therefore we can 
insert the summations "^Xi ' ' ' Sajv+i ^^^^ matrix. Now we arrive at the following. 

Proposition 4.1 The vacuum expectation values of two fields have the determinant for- 
mulas 

,(^Af(/^l, . . . ,/XAr)|V'(xi,0)V'"l"(3;2,t)|^Ar(;Ul, . . .,fiN))e 
e(^'Af(/Ul, . . . ,IJ,n)\'^n{IJ'1, ■ ■ ■,IJ'N))e 



^-ith 



^ ^ ee^^^'c,[x,\s)C,{x2\s) + ^ 



- \^ ^ sGfN 



Dynamical Correlation Functions for an Impenetrable Bose Gas 



111 



e'*'^Ce{x2\s)Je{xi\s,fXk) 



a=0 



Here we have used 



and functions Ce{x\s) and Ie{x\s, fi) are defined in (3.^) and respectively. 

The size of the above matrix depends on the state number N, however, the element of 
the matrix does not depend on A^. By calculations, we obtain 



^ ee^'^'c,{x,\s)C,{x2\s) 



}_ \ ^ its'^—is{xi—X2) 



2L 



I \ ^ its'^-is(xi+X2) 

^ 2L ^ 



(4.1) 



2L 



Y e^''"C,{x^\s)J,{x2\s,^,) = ^ 



2 /vr 



Ats"^ —isx\ 



S — n 



■ sin(x2(s - li)) - (^^-'^^ U e{/i ^ (-M)} 



(4.2) 



and 



2 /vr 



-J e^*^ J,(xi|s, /.)J,(X2|., A) = d,,^ --(-)-= 



1 

\ — [i 
2 /vr 



e**^' sin(xi(A - /i)) + e**^^ sin(a;2(A - ^)) 



(4.3) 



TT VL 

1 



+e- 



E 

^ (-m)} 



s — A s — /Li 



sin((s — ^)xi) sin((s — X)X2) 



X + iJ, 

It is straightforwards to take the thermodynamic limit of the right hand side of the equa- 



tions (|4.1|) , ( [4. 2D and ([4.3|). We arrive at Theorem 2.2. Next we consider the finite tem- 
perature thermodynamics. By statistical mechanics arguments, at temperature T > 0, the 
thermodynamic equilibrium distribution of the spectral parameters is given by the Fermi 
weight T?(A) (PD: 



lim I — 



1 



LJ A7+1 — Xj 
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Therefore the density is given by 



oo 



DiT) = j = l I mdX. 







Now we arrive at Theorem 2.3. 



5 Differential equations 

In this section we will study the most interesting case {tlj{0,0)tlj'^ {x,t))e,T, which only ap- 
pears in open boundary model. For Dirichlet boundary case e = — , (^'(0, 0)^^(x, = 0, 
because the wave functions become zero. We will consider Neumann boundary case 
(^'(0, 0)^^(x, i))+^T and derive the differentaial equations which describe the dynamical 
correlation functions. 

5.1 Preparations 

First we will consider the zero temperature and general xi,X2 case. 
By using the relation: 



d / 2 ~ 
— det 1 Ve - aA, 

Oa \ TT 



det ( 1 - -v}l Ti( (l- -v}l A 



a=0 



vr / \ v vr 



we obtain the following formulae 

(V'(xi,0)V'^(x2,t))e,o = e-'^^^det (^1 - 



-1 



X I G{xi - X2) + eG{xi + X2) - ^Tr ( M - ^V, ] A, 



Define the integral operator by 
q 



Rj) (A) = i?,(A,/i)/(/x)d/x. 



The kernel function i?e(A, /i) is characterized by the following integral equation: 



1 - -VA ( 1 + -Re]= 1. 

vr / V vr 



Define the integral operator S and L by 

(sf) (A) = r s{x, is)f{is)dt,, (If) (A) = r L(A, ^i)f{^i)d^,, 



where kernel function L{X, fi) is defined in ( |1.1D . The kernel function S{X,n) is character- 
ized by the following integral equation: 



1--L) (1 + -5) =1. 
vr / V vr 
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Lemma 5.1 The kernel functions are related by the following linear relation 
R,{X,fi) = S{X,fi) + eS{X,-fi). 



Proof. The following characteristic relation holds: 

5(A, /i) - - r (L(A, u)S{u, fi) + L(A, -iy)S{-u, fi)) du = L{X, ^jl). 
71" Jo 

Using the relations = 1, (e = it) and L(A, — /x) = L{—X,fi), we obtain the following 
characteristic relation: 

{S{X, ii) + e5(A, -fi)) - - r (S(A, u) + e5(A, -i/)) {L{iy, fi) + eL(z., -^)) du 

71" Jo 

= L(A,;u)+eL(A,-/x). 

Q.E.D. 

By using lemma ^.1| , we obtain 
Here we have used 

Uf) (A) = r [/(A, fi)fifi)dfi, where C/(A, ^) = P(A|xi, X2)P(/u|x2, (5.1) 

Here we have used 

[A?yf^ (A) = /(-A). 

In the sequel of this section we will consider the special case that xi = 0, X2 = x and 
e = +. The following simplication occurs: 



-^^(A,m)Li=o. 



X2=X 



62 



»f(-A2+M^) Sm(x(A-/i)) 



Therefore det ( 1 V4 

vr 



det I 1 - -V4 
vr 



A — /i 

dose not depend on time variable t: 



det 1 - -W 



IT 



xi=0,X2=x 

where the operator W is defined in ( |2.8| ). By using the relation: 

P(-A|xi,X2) = P(A| - Xi,X2), 
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we obtain the simplication: 

We arrive at formulae 
(V;(0,0)V^(x,0)+,o 

-iht 



2 TV ( \l + -S\U 

TT 



X1=Q,X2=X 



2e— det ( 1 - ( - i-IV ( (l + ^5 ) [/ 



Xl=0,X2=X 



5.2 Differential Equations 

In this section we will find the partial differential equations of variables t and x. By 
discussion in the previous subsection, it is enough to consider the factor: 

G(x,+x.)-i^TY((l + ^^)f/ 

It is convenient to consider the problem in more general situation. We introduce the 
auxiliary functions Gp{X) and the auxiliary vectors (A) and Cpdi) defined by 



1 r°° 1 

Gp(A) = — P.V. / eris\y2p-y2p-l ,t2p-t2p-l)^g^ 

2ir 7_oo s- X 



_^it2p-iX^-iy2p-iX ^it2p-iX'^ -iy2p-iXQ^(^^^ ^ ^ 



2 / e-**2pM'+j2/2pMGj,(;x) 



TT 



Let us set the integral operators Kp by 
(Rpf) (A)= r KpiX,i,)f{,,)df,, 
where we have used the kernel defined by 



TT 1 



-e^iX)e^{l^). 



2 A-/X 
Let us set 

E^{X) = {E^{X) E^{X) E^{X) EiiX)) 



,{{X) ((l + fK,)ef)(A) ), 



V E^i^i) ) 



(«f(l + f^2))(/.) 
ef(M) 
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By direct calculations we obtain the closed differential equations of the above vectors: 



d_ 



Rl 



dQ 



Ijin) = nPj + [Q,Pj], rnjifi) = -filjifi) + 



Here we have used the matrix Pj defined in ( |2.11| ) and 



Q 



2^ 7-00 J- 



\ 



1 

2^ 



— oo 
oo 



ef (s)e2 (s)(is 



where r(s|y,t) is defined in ( |2.5D . Define the integral operator M by 

Mf) (A) = ^V(A,/i)/(/x)a!M, 
where we have used the kernel defined by 

By direct calculations we obtain the following Propositions. 

Proposition 5.2 The kernel L{X, fj,) is the special case of the kernel M(X,fi) 



-Xi Xi -X2 X2 

t t 



L{X, fi\t,xi,X2) = M [ X,fi 



(5.2) 



(5.3) 



Define the vectors F^{X) and F^(/i) by the integral equations 



F^(A)= Ff(A) FiiX) Fi^iX) F^(A) 



Fi'il^) 



l-=-M] E"^] (A), 



E^ [1- -M 



By usual calculation procedure described in [10|, we obtain the closed differential equations 



d 



F«(^) = L,(/x)F^(^), 



F^(/i) = M,{f^)F^i^,), 



Lj{ii)= tiPj + [b,Pj], 



dt, 



db 



(5.4) 



M,(//) = -/xL, (//) + — 
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Here we have used matrix b defined by 

b = B + Q, (5.5) 
where Q is defined in (|5.2| ) and the matrix B is defined by 



The compatibihty condition of the above differential equations (5.4) yields the differential 
equations ( 2.10| ). 



Proposition 5.3 A factor of correlation functions can be written by an element of the 
matrix b 

o(....,,-±^.((,.|5)g)...,(7 ; 7 - 

Proof. The kernel U{X,^) ( [5l| ) is related to the vectors E^{\) and E^[ii) 

Q(A,M) = -2vri^f(A)ii;f(M)|^^_^^_,^^^3_^^_^^,^^,^^o,3=*.=f 
By using Proposition |5.2|, we arrive at the following 



Q.E.D. 



Now we arrive at Theorem |2.4 For finite temperature case T > 0, we prepare some 
functions. Let us set 



(Mt/)(A)=/ M{\^l)^{|x)d^l, 



where ??(//) is defined in (1.2) and M{X,fi) is defined in (^.3|). Define the vectors F^{X)t 
and F^{^)t by the integral equations 

Define the matrix b'^ by 

b^ = B^ + Q, (5.6) 
where Q is defined in ( |5.2| ) and the matrix i?^ is defined by 

«?i=^i'( ?; s s ) =/>''(A).^.^(AMA. 

By the similar discussion as temperature T = case, we arrive at Theorem 
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6 The time-independent Case 



The purpose of this section is to derive the Fredholm minor determinant representations 
for finite-temperature fields correlation functions: 

Our Fredholm minor determinant representations coincide with the one which has been 
obtained in M. When we take the limit t — > 0, the following simplifications occur: 



P{\\xi,X2) e""i^. 
Therefore we obtain 



1 



X — fi 



(sin(xi(A - /i)) + sin(x2(A - fi))) 



:r-(a;i,X2) 



mxMx2))e,T = ^ ) 1 - ^^^^T - aW.^T 



Q = 



1 



2^1 



.-det(l--y. 



vr 



Tr 



1 - -V,,T 
vr 



Here the integral operators are given by 



Ve,Tf]{\) 



K,r(A,^)/(/i)d/i, 

where the integral kernels are given by 

1 



W^.,T / (A) 



A — /i 



{sin(2;i(A - /i)) + sin(x2(A - /i))} 



+ e- 



{sin(xi(A + /x)) + sin(x2(A + /x))} 



Pay attention to the Fourier transforms: 



/(A) 



1 



The following identity holds for functions /e(eA) = e/(eA): 

) 



dXy^e-'^^fiX). 



— sin(x(A - /i)) + e— — sin(x(A + //)) 
X — fi A + /i 



where 



'd{v) {cos((^ — rfjv) + e cos((^ + r])i^)} du. 
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Therefore we arrive at 



det I 1 - -V^T 

IT 



det 1 



vr 



where the integral operator 6^^'^^^ is defined by 

[eiy^f) io = I ((E(yi - + E(y2 - aweH^, a) fiOdi'. 

Here E(^) represents the step function 



1, for ^ > 0, 
0, for ^ < 0. 



Let us set 



Ae(e,r?) 



det 




e J 


det 






) 





Here we have used the fohowing notation of the r-th Fredholm minor determinants: 
m ■■■ r]r 



det 1 - XKi 



+r 



n 



n=0 

where we have used 

6 ••• e 



dX 



1 • • • 



dXnK, 



n+r 



6 



TJr Xl 



Xn 
Xn 



Kr, 



m 



m 



, det 



The integral operator Kj is defined by using the integral kernel K{X, fi) and the integral 
interval I: 

ikif){X) = J^K{X,f,)fif,)dii. 

From the above definition of the Fredholm minor determinants, the function Ae(^,?7) 
satisfies the integral equation: 

J -XI 71" 

Let us take the Fourier transforms of this integral equation 



1 



Jo 



27r0%) 



TT 

Therefore we obtain 
A<(?,>() 



vr 



1 - -Ve.T 
vr 



W. 



Now we arrive at Theorem p. 61. In much the same way as with finite temperature case 



T > we arrive at Theorem 2.7. 
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